Strong gate coupling of high- Q nanomechanical resonators 
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The detection of mechanical vibrations near the quantum limit is a formidable challenge since the 
displacement becomes vanishingly small when the number of phonon quanta tends towards zero. 
An interesting setup for on-chip nanomechanical resonators is that of coupling them to electrical 
microwave cavities for detection and manipulation. Here we show how to achieve a large cavity cou- 
pling energy of up to (2tt) 1 MHz/nm for metallic beam resonators at tens of MHz. We used focused 
ion beam (FIB) cutting to produce uniform slits down to 10 nm, separating patterned resonators 
from their gate electrodes, in suspended aluminum films. We measured the thermomechanical 
vibrations down to a temperature of 25 mK, and we obtained a low number of about twenty phonons 
at the equilibrium bath temperature. The mechanical properties of Al were excellent after FIB 
cutting and we recorded a quality factor of Q ~ 3 x 10 5 for a 67 MHz resonator at a temperature of 
25 mK. Between 0.2K and 2K we find that the dissipation is linearly proportional to the temperature. 
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The measurement of small-amplitude vibrations in me- 
chanical systems is becoming an increasingly interesting 
problem [1, 2]. From the point of view of basic science, 
the study of mechanical systems close to the quantum 
limit has attracted a lot of interest recently [3, 4]. The en- 
deavor towards the ground state of the harmonic phonon 
oscillations has been going on in various physical systems 
such as in optomechanics [5-7], or in electrically cou- 
pled beam resonators which have been measured using 
single-electron transistors [4, 8], or lately, with on-chip 
microwave cavities [9-12]. 

The quantum challenge is posed by several issues, in- 
cluding the relatively low frequency (/o ~ 10 MHz), 
of the lowest modes in suspended beams. The quan- 
tum limit implies stringent requirements on temperature, 
since hfo needs to be small in comparison to ksT . On the 
other hand, at higher frequencies, the coupling to mea- 
suring systems diminishes rapidly. Third, the zero-point 
vibration amplitudes xzp = \/K/2mu)Q, where m is the 
effective mass and ujq = 2wfo is the angular frequency, 
are vanishingly small even at the atomic scale. Very re- 
cently, O'Connell et al. [13] demonstrated a piezoelectric 
mechanical mode at the quantum ground state by using a 
coupling to a superconducting qubit. However, bringing 
a purely mechanical mode to the quantum limit remains 
an ongoing quest, with the goal becoming a reality prob- 
ably in the near future. 

Micromcchanical resonators are also used in applica- 
tions as detectors. The best devices take advantage of the 
trend to smaller size and higher frequencies, and will soon 
approach sensing at the atomic mass unit level [14, 15]. 
They could also be operated as sensors of position, force, 
or high-frequency electromagnetic fields. 

For conductive resonators, capacitive coupling to an 



electrical measuring apparatus is useful for readout. In 
contrast to magnetomotive or optical detection, one can 
obtain very low electrical losses. An impediment to ca- 
pacitive measurement is the reduction of coupling with 
increasing frequency that accompanies reduction in size. 
The coupling via a capacitance C g is proportional to 
dC'g/dx, which for typical geometries depends on the vac- 
uum gap xq between the vibrating beam and the gate 
electrode as (xq)~ 2 , and hence small gaps are crucial. 
Using traditional e-beam fabrication, it is hard to ob- 
tain Xq below 100 nm. Here we show that focused ion 
beam (FIB) milling yields gaps down to ~ 10 nm for res- 
onator beams up to 10 micron in length and 100 — 200 nm 
in thickness. An improvement in coupling by two orders 
of magnitude can thus be achieved using FIB-tcchniqucs. 
We also find that mechanical dissipation is reduced by the 
FIB milling, and we observe quality factors approaching 
half a million at dilution refrigerator temperatures. 

Our system, as shown in Fig. 2, is a nexural beam 
coupled capacitively to an LC cavity resonator detuned 
substantially from the beam's frequency, i.e., /lc 3> fo- 
The beam displacement x pulls the cavity frequency 
by the amount gx, where the coupling energy is g = 
(cj lc /2C) (dCg/dx), and w LC = 27r/ LC . In the quan- 
tum language, the Hamiltonian of the whole system is 
H = hw^Qcva + Huj(,b'b + Hgxa/a — igv(b + b )(o — 6 )■ 
Here, d^ and a are the creation and annihilation opera- 
tors of the cavity, and similarly w and b for the nanores- 
onator. x = xzp + bj is the quantum amplitude of 
the nanoresonator. The photon-exchange linear coupling 
with the energy gy — xzpV g \J~^f- (dC g /dx) can be 
tuned up to a very high value with a dc voltage Vdc, 
but is negligible unless the cavity and the nanoresonator 
are close to resonance. 
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FIG. 1: The fabrication of strongly capacitive coupled nanoresonators, (l)-(4), on top of a high-resistivity silicon wafer. E-beam 
lithography was used to pattern the metallization. Before Focused Ion Beam (FIB) milling, the e-beam evaporated aluminum 
T-structure was suspended by an optimized SF6 Reactive Ion Etch. The single-pass line cut, with current 10 pA, milling depth 
~ 100 nm, and ~ 5 nm at 30 kV, gives vacuum gaps down to 10 nm without any noticeable bending of the beam. On right, 
SEM image of a beam-gate structure (sample C). The detail shows a short beam with a < 10 nm gap, and typical rounded 
cross section seen after milling. 



We have used [11] shorter beams, I ~ 1.5. ..10 fim, than 
in Refs. 10, 12 and we obtain higher frequencies [16]. 
Therefore, one would obtain low mechanical occupancies 
already at dilution refrigerator temperatures. The trade- 
off, however, is the smaller area for detection due to the 
shorter beam length when measuring via an LC cavity, as 
well as a reduction of gain as oc /lc//o- To this end, FIB 
processing offers a well suited method to make extremely 
narrow vacuum gaps for relatively short beams. 

In our fabrication, we used focused ion beam (FIB) 
milling to separate the gate and the beam from a prefab- 
ricated, T-shaped, suspended structure, Fig. 1. The pre- 
fabrication followed closely the steps presented in Ref. 11. 
A brief anneal at 250 C for 30 min was performed to 
generate a built-in tensile strain before cutting. In order 
to avoid contamination, the FIB modification was per- 
formed without exposing nearby structures to the Ga- 
beam. This was achieved by dual beam operation with 
careful across-beam alignment (Dual Beam Microscope: 
Helios Nanolab, FEI Ltd.). The line cut was performed 
at 30 keV in single-pass mode, instead of gradually scrap- 
ing the surface, in order to exploit differential sputtering 
effects which allow smaller gap dimensions and signifi- 
cantly reduce Al redeposition onto the structure. Vac- 
uum gaps of 50-10 nm were obtained by an appropriate 
choice of the ion beam diameter (min 0^5 nm). The 
dimensions of the measured samples are listed in Table I. 

An interesting issue in the fabrication process is that 
of possible side-effects of FIB milling on the sidewalls of 
the slit. These include Ga implantation, amorphization, 
heating, and generation of defects. Investigations of these 



TABLE I: Length I, width t x , height t y , and the vacuum 
gap xo of the measured nanomechanical resonators. The res- 
onance frequencies /o were measured at 4.2 K in a vacuum 
better than ~ 10 -4 mBar. The estimated frequency is given 

by fo = 7*\/fV /l + £ T^T'' where E > 1 = ^i!' L > and P are > 
respectively, the Young's modulus, the second moment of in- 
ertia, and the density, e = Al/l ~ 0.6% denotes the tensile 
strain. /o h and /q XP denote the calculated and measured me- 
chanical resonance frequencies, respectively. 



sample 


i(/im) 


t x (nm) 


£ H (nm) 


io(nm) 


/,5 h (MHz) 


/o Cxp (MHz) 


A 


4.2 


125 


141 


24 


65 


67.3 


B 


3.0 


147 


141 


29 


114 


101.3 


C 


5.4 


161 


143 


18 


46 


50.8 


D 


5.3 


187 


143 


21 


52 


52.9 


E 


9.3 


309 


155 


20 


27 


27.2 


F 


9.3 


290 


155 


14 


27 


26.7 



artifacts have mainly concentrated on substrate materi- 
als and semiconductors [17, 18]. However, some general 
conclusions can be made on the basis of our process pa- 
rameters. The thickness of the sidcwall layer with dam- 
age is on the order of 10-30 nm for the employed beam 
energy of 30 kV [17-19]. Since this is not negligible in 
comparison to the width of our beams, both mechani- 
cal as well as electrical properties may be altered by the 
FIB-induced damage. 

The damage is known to lead to stress enhancement in 
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FIG. 2: (a) Reflection measurement scheme for the detection of nanomechanical vibrations coupled capacitively to the electrical 
resonance frequency of a cavity (lumped element LC circuit in the figure), /lc = fLc(l + gx)- The vibrations of the beam at the 
mechanical frequency fo <C /lc — (£C) _1//2 /27r change the total capacitance. The resulting sideband voltage is measured by a 
spectrum analyzer. The inductive coupling of the cavity allows for introducing a dc voltage to the gate; (b) electrically actuated 
response of a 5.4-/im-long nanomechanical resonator (sample D) as a function of mechanical drive frequency f g measured at 
T — 4.2 K, with voltages Vdc = 0.2 V and Vlc — 1 niV, while the drive for the mechanical mode Vi is increased from bottom 
(60 nV) to top (6 mV) at 10 dB intervals; (c) resonance measured as a function of dc bias Vdc (sample C). The solid curves 
display the calculated dependence at gap values xq indicated in the figure. 



the milled structures, due to either formation of solid so- 
lution, or Taylor or precipitation hardening [19—21]. Ion 
doses required to cut the film led to bent beams, an in- 
dication of compressive strain induced in the surface ex- 
posed to milling. Stress effects were negligible for beams 
of relatively low aspect ratio l/t x < 30, confirmed by the 
uniform structure in Fig. 1. We use differential thermal 
contraction to estimate the additional stress on cooling, 
T - Et x t v ^ with - 4 • 1CT 3 . Together with E = 70 
GPa and an initial stress of 0.15 GPa, we obtain resonant 
frequencies within 10% of those listed in Table I. 

We first detected and characterized the flexural vibra- 
tions by electric actuation through the gate capacitance, 
see Ref. [11] as in Fig. 2. In order to maximize the cou- 
pling, the stray capacitance C in parallel to the beam 
should be minimized. For the best setup (sample E), we 
used an on-chip spiral coil inductor with 15 turns and 
3 jj,m lincwidth. The stray capacitance was further re- 
duced by removal of a substantial part of the underlying 
substrate by the Si etch. The values C ~ 30 fF, L = 14.5 
nH were obtained from electromagnetic simulation of the 
structure, and which agree with the measured /lc = 7.64 
GHz. For the other samples (apart from E), we used off- 
chip inductors with C ~ 0.3 pF, and L ~ 3 — 10 nH. 
The cavity linewidth is adjusted by the tuning capacitor 



C T ~ 1 - 5 pF such that A/ LC < fo- 

The force driving the beam is given as F(t) = 
\Vgtt)~dx~- The mechanical motion gives rise to a time- 
varying capacitance Cn(t) = C g o + f^jfj x. The vari- 
ation of the total capacitance of the LC resonator leads 
to its frequency modulation, which creates sidebands at 
f± = /lc i fo in a microwave reflection measurement. 
The blue sideband voltage V+ was directly measured by 
a spectrum analyzer as shown in Fig. 2(a). As a first 
amplification stage, we employed preamplifiers with a 
noise temperature TV varying between ~ 3. ..300 K. Two 
circulators in cascade were employed to reduce pream- 
plifier backaction noise. The low-temperature measure- 
ments were performed in a dilution refrigerator that also 
provides an ultra high vacuum environment. 

Fig. 2(b) displays the driven response for a 5-/im-long 
Al beam with a gap of 21 nm (sample D in Table I). 
The peaks were measured by sweeping f g and stepping 
Vi by 10 dB between successive spectra. We observe a 
linearly increasing Lorentzian response, which displays a 
Dufhng-type instability (stiffening spring) at deflections 
on the order of 2 nm as expected for a doubly clamped 
beam resonator. 

In Fig. 2(c) we display the measured side band res- 
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FIG. 3: The total energy of the thermomechanical vibration 
peak measured for sample E (/o = 27 MHz) as a function of 
the cryostat temperature. Linear behavior is observed down 
to the base temperature of 26 mK, where we obtain a low 
number of quanta n ~ 20 without electrical cooling via the 
cavity. The insets show the thermal vibration spectra mea- 
sured at two different temperatures. The red lines are fits 
(see text). The width of the peaks changes due to the strong 
temperature dependence of the mechanical quality factor. 



onance position as a function of dc bias Vd c . The 
concave parabolic shape with gate voltage is charac- 
teristic of resonators with small tension. Since the dc 
voltage dependence becomes very sensitive to the vac- 
uum gap Xq at small separations, this dependence can 
be used for verification of the gap size. Indeed, we 
obtain the best fit for a gap which is close to that 
seen in an SEM image (18 nm), whereas larger val- 
ues are clearly excluded. The solid curves display the 
calculated dependence modeling the capacitance C g of 
two parallel beams with circular cross section C g = 



tteqI < In 



This cs- 



X /ty + 1 + y (x /ty) +2x /ty 

timate gives, for instance, C g ~ 0.5 fF for sample E. We 
also infer that the FIB processed regions remain metallic. 

Attempts to cool mechanical systems down to the 
quantum ground state have usually taken place at rel- 
atively low frequencies, using active feedback cooling, or 
by passive sideband cooling [5, 10, 12]. Using the latter 
method, occupation quanta number ~ 4 were recently 
obtained [12]. At higher frequencies, electrical cooling 
is not necessary to achieve low quantum occupation, as 
demonstrated by the Cleland group [13]. For our res- 
onators around 50 MHz, relatively low occupancies are 
feasible already without sideband cooling. 

We measured the thermomechanical vibrations at sev- 
eral different temperatures for the 9 /im long beam, sam- 
ple E. The temperature was simultaneously monitored 
by a primary Coulomb blockade thermometer. The noise 
floor is set by the cold amplifier noise, Tjv ~ 2.5 K. 



We fitted a Lorentzian (representing the thermal vibra- 
tion peak), plus the noise floor to the measured sideband 
power, see insets of Fig 3. The area A under the peak, 
given by the parameters of the fitted Lorentzian, is di- 
rectly proportional to the energy of the vibration degree 
of freedom, but requires knowledge of the system gain 
which is difficult to deduce accurately. At higher tem- 
peratures, we assume that the beam is thermalized and 
calibrate the area of the peak against fc^T = A = nHcoo. 
We obtained linear behavior, see Fig 3, of A against T 
down to the base temperature of 26 mK, and hence infer 
a low number of quanta ~ 20 at 26 mK. 

The coupling energy g between the resonator and cav- 
ity can be estimated using the parameters of the cavity, 
and dCg/dx ~ 10 nF/m based on the dimensions of the 
beam and the gap. The model is verified by the dc volt- 
age dependence, see Fig. 2. We obtain an order of magni- 
tude larger coupling g ~ 2?r X 1 MHz / nm than previously 
demonstrated [10, 12, 22]. 

A high mechanical quality factor Q m is instrumental 
for various detector applications, as well as side band 
cooling near the quantum ground state, and to observe 
phenomena in that regime. The losses are often dom- 
inated by defects which influence the mechanical and 
thermodynamic properties of amorphous materials at low 
temperatures [23], described using two level tunneling 
states (TLS). Typically, the quality factor decreases with 
the physical size of the resonator, an indication that sur- 
face losses pose a bottleneck [2]. A suitable figure of 
merit to take into account both frequency and Q value is 
their product, foQ m - For sample A, we obtain, with the 
lowest mode at 67 MHz and Q m = 2.7 x 10 5 at T = 25 
mK, f Q m ~ 2 x 10 13 Hz. 

Interestingly, we observed that the Q value seemed to 
be improved by the FIB milling. A test sample fab- 
ricated using a single e-beam lithography step, having 
I = 2.5 /im, fo — 105 MHz, displayed an order of mag- 
nitude lower Q m = 25 x 10 3 at T = 25 mK. The im- 
provement of Q can be due to modification of the proper- 
ties of microscopic defects. A more complicated scenario 
could be envisioned owing to the compressive strain gen- 
erated by the FIB cutting at the edges of the beam. This 
would cause the edges to have less strain than the rest 
of the beam. Therefore the stress due to flcxural motion, 
which is usually concentrated near the beam edge would 
be more uniformly distributed over the beam. The cou- 
pling of phonons to defects depends on the strength of 
the local strain field, and thus is reduced. 

Wc also investigated the temperature dependence of 
the performance of some of the samples. We used the 
driven response, as in Fig. 2(b), or a different approach 
which uses two microwave frequencies /lc and ji such 
that /a — /lc ~ fo [26]. Both methods gave similar re- 
sults. In the latter case, the low-frequency force driving 
the mechanics is now oc Vlc^2i where VLc an d Vi are the 
microwave voltages, fi was kept constant, while rcflcc- 
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FIG. 4: Measured temperature dependence of the inverse me- 
chanical quality factor (samples A and E) . The solid line indi- 
cates oc T dependence. The inset is an example of the driven 
response (sample E), at T = 25 mK, fitted to a Lorenzian 
which yields Q m — 3.64 x 10 s . The right frame shows the 
corresponding change in the mechanical frequency, together 
with a logarithmic fits. 



tion at /lc was swept with a vector network analyzer. 
We find a strong decrease in the dissipation Q m x with 
decreasing temperature in the range 0.1 - 2 K, see Fig. 4. 
Our data follow approximately Q" 1 oc T in this range. 
This behavior is qualitatively similar to that observed in 
latest experiments measured in a considerable magnetic 
field. Much weaker T-dependence was obtained in amor- 
phous gold resonators [24]. Hoehne et al. [25] obtained 
recently a similar linear temperature dependence in litho- 
graphically patterned Al resonators. Therefore, we con- 
clude that while the dissipation is decreased due to FIB 
processing, its source likely stays unchanged. Notice that 
the saturation below 0.1 K is not due to thermal decou- 
pling, since we observe a consistent linear dependence of 
the thermal vibration energy down to 26 mK, see Fig. 3. 

The mechanical resonance frequency /o was also found 
to depend on T, although rather weakly: at V dc = 0.2 V, 
Sf/f ~ 2 ■ 10~ 5 In (T/Tq) while at V dc = 0.8 V, wc find 
Sf/f Q ~ 4 ■ 10" 5 In (T/T ). Since grows with P dc , 
part of the frequency shift can be attributed to changes 
in the dielectric constant of the substrate material. 

The logarithmic temperature dependence of 6f/fo and 
the rolloff in 1/Q m observed here (Fig. 4) are qualita- 
tively consistent to those predicted by TLS and explained 
by the tunneling model [23]. In this model, the former 
is due to the resonant interaction between drive phonons 
and TLS. The rolloff in 1/Q m is observed at the temper- 
ature Tqo ~ 1 K at which the relaxation rate of the TLS 
(t -1 ) equals the mechanical frequency, luqt = 1. This is 
consistent with that observed in Ref. 27 since Tco varies 
as w 1 / 3 . 

The temperature dependence of 1/ Q m for T < Tqo de- 
pends on the dominant relaxation mechanism (electron- 
driven or phonon-driven) and size effects. Electron- 
driven relaxation produces a linear temperature depen- 
dence of 1/Qm for T < Tco- However, here, in the su- 



perconducting state, one expects that phonon-driven re- 
laxation should dominate [28] . Yet the T 3 dependence of 
1/Qm expected for T <C Tco due to relaxation driven by 
3D phonon modes [29] is not developed, perhaps because 
other dissipation mechanisms become important at the 
lowest temperatures. 

The dominant thermal phonon wavelength in Al, ~ 37 
nmK/T becomes comparable to the widths and thick- 
nesses of the beams at 100-300 mK [31]. In this limit, 
the dominant modes are no longer 3D phonon modes [32] 
and can contribute a linear temperature dependence 
[25, 30]. This mechanism thus could dominate relax- 
ational damping in our high-frequency resonators. At 
the lowest temperatures, the resonant contribution to 
1/Qm oc tanh (hu/2kBT) may contribute to the mea- 
sured saturation. 

The large capacitive coupling we have achieved be- 
tween a nanoresonator and an LC circuit can be useful 
for various sensitive detectors. The responsivity of the 
signal voltage to a mass change is dV/dm ~ VQ m /m. 
Four our sample (E) set-up at T = 4 K as in Fig. 2, 
we obtain s m ~ 10" 26 kg/ Viz - 10 AMU/Vlz, which 
would improve the state of the art [14] of mass detection 
by an order of magnitude. For studies near the quantum 
limit, the strong coupling achieved via FIB milling is a 
critical asset since it enables Q-enhanced high-frequency 
resonators. Indeed, for short beams 1 — 2 /im in length, 
we have obtained vacuum gaps approaching 5 nm. These 
findings substantially improve, for example, the efficiency 
of electrical cooling via the LC circuit which is inversely 
proportional to the fourth power of the gap. 
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